The on-top pair-correlation density in the homogeneous electron liquid 
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The ladder theory, in which the Bethe-Goldstone equation for the effective potential between two 
scattering particles plays a central role, is well known for its satisfactory description of the short- 
range correlations in the homogeneous electron liquid. By solving exactly the Bethe-Goldstone 
equation in the limit of large transfer momentum between two scattering particles, we obtain ac- 
curate results for the on-top pair-correlation density g(0), in both three dimensions and two di- 
mensions. Furthermore, we prove, in general, the ladder theory satisfies the cusp condition for the 
pair-correlation density g(r) at zero distance r = 0. 
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I. INTRODUCTION 

The pair-correlation density g(r) is one of the key 
concepts in describing the correlation effects, arising 
from Pauli exclusion principle and Coulomb interac- 
tion, in the homogeneous electron liquid (or gas).Q 
It also plays a significant role in the constructions of 
the exchange-correlation energy density functionals in 
density- functional theory (DFT),|2| since in such con- 
structions the homogeneous electron system is conven- 
tionally taken as a reference system. A great deal of 
theoretical progress has recently been made in giving an 
accurate evaluation of g(r), or the more specific spin- 
resolved pair-correlation densities n aa , (r), with g(r) = 
lbn(r) +mW]-SBIiSllli In particular, g(Q), 
the on-top pair-correlation density, which arises totally 
from <7|x(0) since <7tt(0) = 0, has been well known to play 
a special role in DFT.|lfJ The important implication of 
<?(0) was also realized in many-body theory long ago be- 
cause the random phase approximation (RPA),[lj due to 
its lack of accurate description of the short-range electron 
correlations, yields erroneous negative values for q(0 ) 
when the electron densities are not sufficiently high.jlj 

It is well known that, in many-body theory, the long- 
range correlations can be rather successfully taken into 
account in the RPA, while the short-range correlations 
can be properly described by the ladder theory (LT) . 0, 
111 . ITEl fl^ l In this paper, we attempt to investigate the 
short-range correlations in terms of g(0) in the LT, in 
both three dimensions (3D) and two dimensions (2D). 
In fact, investigations in this direction date back long 
ago , and a great deal of achievement has been made. 
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It is necessary here to give some introduction to the LT. 
The effective interaction V e ff(p, p'; q) in the LT between 
two scattering electrons with respective momenta p and 
p' satisfies the following Bethe-Goldstone equation: [23| 



V e ff(p, p'] q) = v(q) + ^v{q- k) 



(l-n(p + k))(l-n(p'-k)) r ,.x m 
x - Z Veff(P,P ,*■), (1) 



where v{q) is the Fourier transform of the Coulomb po- 
tential, n(p) — 9(kp — p) is the momentum distribution 
in the noninteracting ground state and kp is the Fermi 
momentum, and e p — h 2 p 2 /2m. 

As mentioned above, the RPA gives poor description 
of the short-range correlations of the electrons, especially 
for g(r) as r — > 0. In fact, the results for gj i (r) in the 
RPA violate the following cusp condition: fl4l.ll7l.l24 

MM 



dr 



(d - l)a B 



<m(o), 



(2) 



where d = 3, 2 is the number of spatial dimensions, and 
as is the Bohr radius. It was shown recently [2^ that 
the pair-correlation density obtained from the first order 
perturbation calculation does not satisfy the cusp condi- 
tion either. In this paper, we prove that <?tx( r ) calculated 
from V e f f (p, p'; q) of Eq. satisfies the cusp condition. 
This indicates the reliablity of the LT in the calculations 
of the pair-correlation density at short range. 

The short-range structure of the pair-correlation den- 
sity is determined by the behavior of the effective poten- 
tial V e ff(p, p'; q) at large momentum transfer q. In the 
limiting case, one therefore can approximately replace 
the momenta of the scattering electrons by zero in Eq. 

09. 



v, 



eff 



(0,0; q) =«(«)- 5>(q-k) 



xi^%(0,0;k). (3) 

A frequently used approach to solving Eq. J2J in the 
literature is making the following approximation in the 
Coulomb kernel in the momentum summation: |l?lll8ll2ll 

MM 



t>(q-k) = v(q), q> k 
= v(k), k > q. 



(4) 



With the preceding approximation, an analytical solution 
for V e ff(0, 0; q) was obtained which yields the following 



2 



well-known result for g n (0) in 3D, |15i I29L I3C 

-, 2 



911(0) 



(5) 



where A3 = 2ar s /n with a — (4/97T) 1 / 3 and r s = 
(3/47rn) 1 / 3 /as- A similar result was obtained in 2D.[Ti| 



/o(V4Al 



(6) 



where A2 = r s /v2 with r s = l/y/nriaB in 2D. In Eqs. (0 
and 10, /n(^) is the nth order modified Bessel function. 

In this paper we have managed to solve exactly Eq. 
PJl. i.e., without making the approximation of Eq. J1J. 
Our results for <?|j.(0) are 



in 3D, and 



45(45 + 24A 3 + 4A§) 



2025 + 3105A 3 



1512A§ 



256A 3 



15(64 + 25A 2 + 3A|) 



960 + 1335A 2 + 509Ai + 64A 3 , 



(7) 



(8) 



in 2D. Equations (0) and JBJ) are the main results of this 
paper. 

The paper is organized as follows: In Sect. II, we solve 
Eq. <|3J) exactly both in 3D and 2D. In Sect. Ill, we 
derive analytically the expressions of Eqs. (J7J) and JSJ) 
for (?|x(0). We then compare our results with previous 
ones in the literature in Sect. IV. Sect. V is devoted to 
conclusions. Some technical points on the solutions for 
the coefficients of the large momentum expansions of the 
effective potentials are given in Appendix A. In Appendix 
B, we prove the cusp condition in the LT. 



II. EXACT SOLUTION TO THE 
BETHE-GOLDSTONE INTEGRAL EQUATION 
AT LARGE TRANSFER MOMENTUM 

In this section, we present our solution to Eq. © at 
large momentum transfer q in the effective potential in 
both 3D and 2D. To this end, we denote V e ff(0, 0;q) 
as V e ff(q), and reduce the momenta with unit &f, and 
potentials with v(kp), respectively. We present our solu- 
tion for the 3D case in subsection A, and the 2D case in 
subsection B, separately. 



We expand V e ff(q) in the powers of 1/q, 



2n+2 ' 



(10) 



It can be easily confirmed by iteration that no odd power 
terms in the expansion of V e f / (q) exist in the solution to 
Eq. ©. The erroneous odd power terms introduced 
into V e f f (q) in Refs. 0, |2^, [3(J are purely due to the 
approximation made in the Coulomb kernel in Eq. Q). 
We substitute Eq. <fT~0f> into Eq. ((SJ), and obtain 



EObfi 1 A3 v' — > 



71=0 



q 2 2g 



(11) 



n=0 



where 

M 2n+3 (q) = 



1 „ 1 h 

dk ^Ts 1 ^^ n>0. (12) 



By carrying through partial integration on the right hand 
side of Eq. l(T2*|) . one has, 



1 



2rc + 2 



ln*il-2 9 $ n+1 ( 9 ) 
9-1 



where 



*n+i(g) 



OO j 

dk 



1 



fc 2n+2 k 2 _ q 2 ■ 



(13) 



(14) 



$ n +i(g) defined in the preceding equation can be evalu- 
ated to be 



*n+i(g) 



n 1 

V — 

/ -< n 2m- 



1 



q 2m+2 2(n-m) + l 

m— 



In 



q + i 



2q 2n + 3 ~ q - 1 
Substituting Eq. l(TH)l into Eq. (fTTTfl yields, 



(15) 



2n+3 



(q) = 



1 



n + 1 
1 / 1 



n , 

„2m+l 



q 2m+i 2( n - m ) + 1 



2\q 



,2n+2 



- 1 In 



9-1 



q + i 



, n > 0. (16) 



Finally, substituting Eq. (|16fl into Eq. I|ll|). and com- 
paring the same power orders of 1/q, one obtains the 
following equations for a n : 



3D 



After carrying out the angular integrations in the sum- 
mation of the momentum k, Eq. @ becomes 



1 A r°° 1 



1 i„ | g + fe | 

q - fc 1 



(9) 



and 



a 



1-AaE 



n=0 



2n+V 



2n 



1=0 



ai 



2(1 - n) + 1 



, n > 1. 



(17) 



(18) 
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Equations 1|17|) and (|18|l for a n can be solved exactly 
in principle. In fact, by making the truncation of a n = 
for n > 3, a nearly exact solution can be obtained as 



where 



a 



45(45 + 24A 3 + 4A 3 ! ) 

Da 



and 



where 



Oi 



('2 



15A 3 (45 + 8A 3 ) 



45A 3 (3 + 4A 3 



D 3 = 2025 + 3105A 3 + 1512A? + 256A?. 



(19) 



(20) 



(21) 



(22) 



In Appendix A, we show that the preceding solution 
for ao, which is directly related to (0), as shown in the 
next section, is very close to the exact numerical solution 
to Eqs. HI 7(1 and 1(18(1 . In fact, the large momentum 
behavior of V e ff(q) is dominated by the leading terms in 
the large q expansion of V e ff(q) in Eq. (JTDJ, and hence a 
truncation solution like the preceding one is almost exact. 



B. 2D 

In 2D, we make use of the following expression, 
2?r 



J dre l(q - k) r i, 



(23) 



|q-k 

and rewrite Eq. © as follows: 

Veffiq) = ~q~ttr I dke{k ~ 1] h Veff{k) 

x J dre 4(q - k) r i. (24) 

Carrying out the angular integrations of k and r, we have 
1 r°° f°° 1 

Veff(q) = - - A 2 J dr J dk-V eff (k) 

xJ Q (qr)J (kr), (25) 

where J n (x) is the nth order Bessel function. We expand 
V e ff(q) in the powers of l/q as follows: 



Veff(q) = Ta^+i- 

n=0 q 



(26) 



No even power terms exist in the solution to Eq. ((25(1 . 
Again, the erroneous even power terms |18| appear in 
V e ff(q) due to the approximation made in Eq. @. 
We substitute Eq. lt26*|) into Eq. and obtain 

oo . oc 



N 2n (q) = / dxJ (x) / dk—J (kx/q). (28) 



Carrying out the integration over k in Eq. pg|l. one 
obtains, 



N 2n (q) =E^ 2 ) 



-i (n-1)! 



(n — m)' 



x / dxJ (x/q)J m (x)/x m , n > 1. (29) 
Jo 

The integral on the right hand side of Eq. 1(29(1 can 
be expressed in terms of the hypergeometric function as 
follows, H3 

dxJ {x/q) J m (x)/x m = — ^ T - 

o 2"T(l)r(m+i) 

x fQ, -m+i;l;^V n> 1,(30) 

where T(a) is the Gamma function. Therefore, one has, 

(n-1)! 

m— 1 



iV 2rl (<z) = £ (-2)' 



(n-m)!(2m- 1)!! 



1 



1 . 1 



x Fl -,- m+ -;l;_ I n >l. (31) 

Substituting Eq. JjEJ into Eq. (57^. and comparing the 
same power orders of l/q, one finally gets 



Co 



(32) 



n=0 



and 



oo I 



1=0 m=0 

(2m + l)(2m - 1) . . . (2m - 2n + 3) 



(Z-m)!(2m + l)!! 



, (33) 



for ?i > 1. 

Similarly to the 3D case, Eqs. ((32(1 and 1(33(1 can be 

solved exactly in principle. In fact, a nearly exact so- 
lution can be obtained as follows by the truncation of 
Cn — for n > 3: 



co = 



ci 



and 



where 



C2 



15(64 + 25A 2 + 3A|) 

30A 2 (8 + A 2 ) 
I>a 

45A 2 (1 + A 2 ) 
A, 



n=0 



q q 



D 2 = 960 + 1335A 2 + 509 \\ + 64A^ 



(34) 
(35) 

(36) 
(37) 
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III. RESULTS FOR g n (r) AT SMALL r 

The spin-antiparallel pair -correlation density in the LT 
can be shown to be fijl [lg 



911 



x y e// (p,p';q) e ^ r | 2 , (38) 

where the prime on the summations over p, p' means the 
restrictions < p,p' < hp, and D(p, p'; q) is defined as, 

^(p,p / ;q)= (1 " n(P + ! ))(1 ~r (p/ ~ q)) - ( 39 ) 

Below we present the results for the 3D and 2D cases 
in subsection A and B, respectively. We will reduce r 
with unit 1/kp. 



A. 3D 

Using the approximate solution V e ff(q) for 
V e ff(p, p'; q), one obtains [l7| 



1 - A 3 / dqV e ff(q)jo(qr) 



9ll( r ) = 
Trivially, 

With the expression of Eq. pO(l , one has 
Su(°) = a o- 



l-As/ dqV eff (q) 



(40) 



(41) 



(42) 



Equation (|17|l has been made use of in obtaining the pre- 
ceding result. The expression for ao is given in Eq. I ll 91 1 . 
with which we obtain the final result of Eq. 0. Fur- 
thermore, it is straightforward to show, from Eq. (|4()Jl . 
that at small r, 



9u( r ) = 9u(°) + 2 A 35u(°) r - 



(43) 



B. 2D 



In 2D, one has, 



9u( r ) 



oo 

1-A 2 / dq-V e ff(q)J (qr) 



n q 

Similar derivation to that in the 3D case leads to 

■2 



Su(°) = Cq, 



(44) 



(45) 



or, by the use of Eq. the final result of Eq. (JHJ . 

Furthermore, from Eq. 1441) . one can obtain 



(1.35 
0.3 
(1.25 

(i: 



c 



(1.(15 




3D 






^\This work 




Ovcrhauscr 




Yasuhara 



FIG. 1: On-top pair-correlation density multiplied by r 3 in 
3D. The result of this work fEq. J/11 is compared with Eq. 
@ (Yasuhara's formula in Refs. (I5ll29ll30j| '). and Overhauser 
formula (Eq. (26) in Ref. Q). 



Evidently, the cusp condition of Eq. J5| is satified 
both in Eq. g3J and Eq. 1(35)1. In fact, in Appendix B 
we shall show that, Eq. is satisfied, in general, in the 
full ladder theory. 



IV. COMPARISONS AND DISCUSSIONS 



First of all, at limiting high density, we have, from Eq. 



0, 



g n (0) = 1-2A 3 

= 1 - 0.663r s , 



(47) 



<m0) = <m(o) + 2A2<7u(Qh 



(46) 



in 3D. Equation (|47l) is the same as the corresponding 
Yasuhara's result. |23] We note that, the first order per- 
turbation calculation, nam nam which is believed to 
approach to the exact result at high density limit, yields a 
result of 1 — 0.73177V We plot g(0) x r s calculated from 
Eq. 10 in Fi g. 1, in comparison with that calculated 
from Eq. JSJ [ljj |2|| |30j . Notice that the discrepancy 
between Eq. and Eq. 0, which appears not minor, 
arises purely from the approximation of Eq. made in 
obtaining Eq. in Yasuhara's theory. In effect, Lowy 
and Brown |16| had thrown doubt on the validity of the 
approximation of Eq. 0. We hence justify their doubt, 
at least for the limiting short range correlations. The 
result based on Overhauser's proposal (Eq. (26) in Ref. 
|33j) is also shown in Fig. 1. The comparison hence in- 
dicates that the coincidence between Overhauser's result 
(and the corresponding numerical result of Gori-Giorgi 
and Perdew |j|) and Yasuhara's is accidental. 

In Fig. 2, we plot <?(0) x r s in 2D calculated from Eq. 
0, together with that from Eq. ©.18] Once again, 
we emphasize that the discrepancy is totally due to the 
approximation of Eq. made in obtaining Eq. 0. 
However, at limiting high density, both equations yield 



5 



c 




FIG. 2: On-top pair-correlation density multiplied by r 3 in 
2D. The result of this work (Eq. JHJl) is compared with Eq. 
(Nagano et al. formula in Ref. Hg|)> an< A Polini et al. 
formula (Eq. (17) in Ref. 0). 



o 0.6 



Exact result 
Eq. (19) 
Eq. (Al) 



FIG. 3: ao calculated from Eq. I19H and Eq. JAIL together 
with the exact numerical result. The inset shows the large r 3 



region. 



the following same result: 

<7 U (0) = 1-2A 2 . 



(48) 



For a comparison, we have also shown in Fig. 2 the re- 
sult of Eq. (17) in Ref. 0, which was proposed by Polini 
et al. based on an interpolation between the first-order 
(second-order in terms of the correlation energy) calcu- 
lation for the weak-coupling limit and Overhauser type 
calculation [33| for the strong-coupling limit. 



V. CONCLUSIONS 

The proper approach to the short-range electron cor- 
relations in many-body theory is the ladder theory in 
which the effective potential between two scattering par- 
ticles satisfies the Bethe-Goldstone equation of Eq. 
In this paper, we have proved that, the ladder theory sat- 
isfies the cusp condition for the pair-correlation density 
in the homogeneous electron liquid. This enhances our 
belief in the capability of the ladder theory in describing 
the short-range correlations, especially in calculating the 
pair-correlation density. 

The main results obtained in this paper are, in effect, 
Eq. and Eq. JSJ) given in the Introduction, in three 
dimensions and two dimensions respectively, for the on- 
top pair-correlation density in the homogeneous electron 
liquid. These results have been derived by solving Eq. 

in which the two scattering particles in the Bethe- 
Goldstone equation are approximately taken to be static. 
This approximation should be reasonable since the limit- 
ing short range structure of the pair-correlation is deter- 
mined by the large transfer momentum behavior of the 
effective potential. The major theoretical progress made 
in this paper is that we have removed the approxima- 
tion of Eq. (0J frequently made in the literature for the 



Coulomb kernel in solving Eq. Our solution to Eq. 
(0 is thus exact. 
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APPENDIX A: ON THE SOLUTIONS TO EQS. 
ifTTjl AND (HHJl, AND EQS. (f32j> AND (|33|). 

A nearly exact solution for ao to Eqs. (|17|l and Ijl8(l 
has been given in Eq. I|19fl in Sect. Ill by the truncation 
of a n — for n > 3. Below we give the solution for a n 
by the truncation of a n = for n > 4. 



a 



175(14175 + 9585A 3 + 2520A^ + 256X 3 3 ) 



(Al) 



105A 3 (7875 + 2480A 3 + 256A§) 
a 1 - ~- , (A2) 



105A 3 (1575 + 2280A 3 + 256A§) ..„. 
«2 = ~ , (A3) 



and 



where 



D, 



175A 3 (405 + 576A 3 + 256A|) . . AS 

a 3 = ^ ~ ^, (A4) 



D 3 = 2480625 + 4158000A 3 + 2437200A 3 

+634880A| i + 65536A|. (A5) 
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FIG. 4: co calculated from Eq. 13411 and Eq. 1A6L together 
with exact numerical result. The inset shows the large r a 
region. 



In Fig. 3, we plot the results for ao calculated from Eqs. 
(|TU|l and (|A1|) , together with the corresponding exact nu- 
merical solution to Eqs (|17fl . and l|18fl . There is basically 
no difference among them. We present the expressions 
of Eqs. (1a"2|) . (|X3|) . (|A^|) above for possible future 

reference. 

Similar expressions for the 2D case are given below. 



Co = 



35(12288 + 6000A 2 + 1121A| + 80Xl) 



D-2 



Cl 



C2 



and 



where 



c 3 



1680A 2 (64 + 14A 2 + A 2 ) 

w 2 ■ 



105A 2 (192 + 207A 2 + 16A^) 



350A 2 (24 + 25A 2 +8A^) 

w, ■ 



(A6) 



(A7) 



(A8) 



(A9) 



D 2 = 430080 + 640080A 2 + 290307A 2 ! 
+55472A^ + 4096At 

The corresponding illustration is given in Fig. 4. 



APPENDIX B: CUSP CONDITION IN THE 
LADDER THEORY 

Due to the singularity of the Coulomb potential be- 
tween electrons, the many-body schrodinger wavefunc- 
tion has a cusp when any two electrons coalesce. [3 HM 



tThis fact leads to the cusp condition of Eq. (0) 
[Hill El |H H3 for the pair-correlation density 
(also known as Kimball relation in the literature of many- 
electron theory). Recently it was claimed that Eq. @ 
is not satisfied in the LT.|37j In this appendix, we give a 
rigorous proof for Eq. J5J in the LT. The proof will be 
formulated in 3D. 

We start with the definition of the spin-parallel static 
structure factor as follows: 

Sub) = ^t(-q)%(q)> - y<W)> (bi) 

where 7v(q) is the spin-resolved density operator and 
N is the particle number. It has been shown that the 
spin-antiparallel static structure factor in the LT can be 
expressed in terms of the effective potential of Eq. as 



2D(p,p';q)F e// (p,p';q) 
|^(p,p';k)%(p,p' ; k) 

k 

x D(p, p'; k - q)V e ff (p, p'; k - q) 



.(B2) 



Next we examine the large momentum structure of 
S^i(q). For p,p' < kp and q — > co, one has, from Eqs. 
J3J and 



£>(p,p';q) V e ff(p,p';q) = -—v(q) 



1+ ^fl(p,p';k)%(p,p';k) 



(B3) 



which evidently goes to zero in the order of 0(l/q 4 ). 
Therefore 



£ £>(p,p';k)V e// (p,p';k) 



xO(p,p';q-k)y e// (p,p';q-k) 
1 



-v{q) 



1 + E£>(p,p';k)V r e// (p,p';k) 



£ J D(p,p';k')K//(p,p';k') 



(B4) 



In obtaining Eq. I|B4(1 . we have used the following rela- 
tion, 

lim V/(k)/(q-k)=2/(q)V/(k), (B5) 



if limg^oo /(q) ~ 0(1/ q 4 ). It seems that a mistake oc- 
curs in Ref. [23] due to a possible miss of the factor 2 on 
the right hand side of the preceding equation, as it was 
employed to derive Eq. (12) from Eq. (11) in Ref. [3^| . 
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Substituting Eqs. (|B3f) and (|B4(I into Eq. (|B2|I . one 

has 

sn(*) = --^£'[i+ E^p'; k ) 



P.P' 



xK//(p,p';k) 



(B6) 



On the other hand, from Eq. 1)38(1. we have 
.9u(0) = 4E'[ 1 + E ^P'P'JI) 



p,p' 



xK//(p,p';q) 



Comparing Eq. I)B6(I and Eq. 1|B7|) yields 

27re 2 nm 
*%WJ = ^2 ffUl )' 



(B7) 



(B8) 



Combining the preceding result with the following well- 
known relation \2M. 



lim q 4 S n (q)=-2nn^^-\ r=0 , 

q—>oc ijr 



(B9) 



one proves the cusp condition of Eq. (0 in the LT. 

The above proof can be straightforwardly extended to 
the 2D case. In fact, in 2D, it can be similarly shown 



: <m(o), 



(BIO) 



in the LT. Combining the above result with the following 
relation [38| 



lim q S u (q) = --irn — 

q^oo 2 or 

leads to Eq. (J2)l for the 2D case. 



r=0, 



(Bll) 
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